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Abstract 

We present a number of arguments to demonstrate that a quantum analog of Cherenkov 
effect occurs when two dispersive objects are in relative motion. Specifically we show 
that two semi-infinite plates experience friction beyond a threshold velocity which, in 
their center-of-mass frame, is the phase speed of light within their medium. The loss in 



t^ mechanical energy is radiated away through the plates before getting fully absorbed in 

^ , the form of heat. By deriving various correlation functions inside and outside the two 

'— ^ plates, we explicitly compute the radiation, and discuss its dependence on the reference 

1-H frame. 

> 

o 

On 1 Introduction 

^, 

^t An intriguing manifestation of quantum theory in macroscopic bodies is the non-contact 

^_J friction between objects in relative motion. For example, two surfaces (or semi-infinite 

J> plates) moving in parallel experience a frictional force if the objects' material is lossy [1-3]. 

S^ The origin of this force is the quantum fluctuations of electromagnetic field, within and 

Cd between the objects; the same fluctuations also give rise to Casimir/van der Waals forces. 

In brief, quantum fluctuations induce currents in each object, which then couple to result 

in the interaction between them. For moving objects, a phase lag between currents leads to 

a frictional force between them. 

For parallel plates, the friction force is related to amplitude of the reflected wave upon 
scattering of an incident wave from each surface (formalized into a reflection matrix below) 
[1,2]. Due to its quantum origin, friction persists even at zero temperature, where it is 
related to the imaginary part of the reflection matrix corresponding to evanescent waves. It 
is usually assumed that the dielectric (or response) function itself has an imaginary part due 



to dissipative properties of the material; this then leads to an imaginary reflection matrix 
and hence friction [19]. However, this is not necessary as, even for a vanishingly small loss, 
evanescent waves lead to an imaginary reflection matrix. 

We consider a constant dielectric function, which fixes the speed of light within the 
medium^. We show that when the velocity of moving plates, in their center- of-mass frame, 
is larger than the phase speed of light in the medium, a frictional force arises between 
them. This is, in fact, a quantum analog of the well-known classical Cherenkov radiation. 
We elaborate on the relation between the friction, and radiation in the gap as well as within 
the plates. 

Quantum Cherenkov radiation was first discovered by Ginzburg and Frank [4] in a 
rather different setup. They argued that when an object (an atom, for example) moves 
inertially and superluminally, i.e. larger than the phase speed of light in a medium, it 
spontaneously emits photons; see Refs. [5, 6] for subsequent reviews by Ginzburg. This 
phenomenon is intimately related to superradiance, first discovered by Zel'dovich [7] in the 
context of rotating objects and black holes: A rotating body amplifies certain incident waves 
even if it is lossy. The underlying physics is that a moving object (atom) can lose energy by 
getting excited. This is because, at superluminal velocities, an excitation in the rest-frame 
of the object corresponds to a loss of energy in the lab frame. Ginzburg and Frank refer to 
this eventuality as the anomalous Doppler effect [4]; see also Ref. [8]. 

Since these unusual observations span several subfields of physics, we find it useful 
to demonstrate the results by a number of different formalisms. We first generalize the 
arguments by Ginzburg and Frank to prove dissipation effects associated with the relative 
motion of two parallel plates. We then use the input-output formalism of quantum optics to 
derive and compute the friction force based on the scattering matrices. An alternative proof 
follows approaches introduced in the context of quantum field theory in curved space-time, 
making use of an inner product to identify the wavefunctions and their (quantum) character. 
Finally we employ the Rytov formalism [9] which is grounded in the fiuctuation-dissipation 
theorem for electrodynamics and well-known to practitioners of non-contact friction. 

To ease computations, however, we consider a scalar field theory as a simpler substitute 
for electromagnetism. The former shares the same conceptual complexity while being more 
tractable analytically. This is particularly useful in expressing complicated Green's func- 
tions with points both inside and outside each plate, or within the gap between them. The 



^The frequency- independence of the dielectric function follows from a vanishing imaginary part, due to 
Kramers- Kronig relations. 



generalization to vector and dyadic electromagnetic expressions should be straightforward 
but laborious. Finally to avoid complications of the full Lorentz transformations, we limit 
ourselves to small velocities — both the relative velocity of the objects and the speed of light 
in their media. 

The remainder of the paper is organized as follows. In Sec. 2, existing formulas are 
used without proof to compute friction, and to discuss the similarities with the classical 
Cherenkov effect. In this section, we elaborate on friction in a specific example. In Sec. 3, we 
consider a general setup, argue for and derive the friction force, as well as emitted radiation, 
in great detail. This section comprises 4 subsections each devoted to one particular method. 
Specifically, we discuss how the radiation within the plates, and in the gap, depend on the 
reference frame. 

2 Friction 

We start with a scalar model that is described by a free field theory in empty space, while 
inside the object a "dielectric" (or, a response) function e is assumed which characterizes 
the object's dispersive properties. The field equation for this model reads 

V2 + e(a;,x)-^j$(a;,x) = 0, (1) 

with e = 1 in the vacuum, and a frequency-dependent constant inside an object. 

We consider the configuration of two parallel semi- infinite dielectric media ('plates') 
in D spatial dimensions, separated by a vacuum gap of size d. For each plate in its rest 
frame, a plane wave of frequency u and wavevector k is refiected with amplitude ('reflection 
matrix') 

" ^ewVc2-k2+ J6^2/c=^-kJ 

where kn is the component of the wavevector parallel to the surface. This result is easily 
obtained by solving the fleld equations inside and outside the plate, and matching the 
reflection amplitude to satisfy the continuity of the fleld and its flrst derivative along the 
boundary. We are particularly interested in friction at zero temperature which is mediated 
solely by evanescent waves. Such waves contribute to friction through the imaginary part 
of the reflection matrices [1-3]. If one plate moves laterally with velocity v along the x axis. 



while the other is at rest, the friction force is given by (introducing the notationdx = da;/27r) 

f=j^ j'^^W^^- |l_e-2|fcxMi?,i?,|2 0(-^ + ^M, (3) 

where Q is the Heaviside step function, k_\_ = . oj'^ /c^ — kn, and L^^^ is the area. Note 
that the reflection matrix of the static plate is given by Eq. (2), but that of the moving plate 
is obtained after Lorentz transforming to the lab frame. We leave the derivation of Eq. (3) 
to the next section, but discuss its implications here. While this equation has been studied 
extensively in the literature, it is usually assumed that the dielectric medium is lossy, with 
a nonzero imaginary part of e. However, even when Ime is vanishingly small, a frictional 
force can be obtained as follows. With Ime ~ 0, the medium can be characterized by the 
modified speed of light vq = cj ^J~k. The only relevant length scale in the problem (aside 
from the overall area L^~^) is the separation d. We can then construct the frictional force 
on purely dimensional grounds as 

where ^ is a function of two dimensionless velocity ratios. Any velocity could have appeared 
as prefactor (with a correspondingly modified function g); we have chosen vq for convenience. 
For small velocities, the dependence on vacuum light velocity c drops out, and 

with g depending only on the ratio of the velocity v to the light speed in the medium vq. 
Interestingly, at small v, only the modified speed within the media is relevant. 
Now note that the Heaviside function in Eq. (3) restricts to frequencies 

(0 <)a; < vkj;. (5) 

Furthermore, the imaginary part of the reflection matrix Ri, given by Eq. (2), is only 
nonzero when uP' jc? — kn < and eiJ^ jc? — kn > 0, which, in turn, imply 

|a;| > -uolkiil > -uolfexl- (6) 

A similar condition holds for the second plate: \ijj'\ > vo\k'^\ with primed values defined in 
the moving reference-frame. For simplicity, we assume that v, vq <^ c, and thus neglect the 
complications of a full Lorentz transformation. Hence, w' ~ cj — vkx and k'^ ^ k^ — Vjjjf? ~ 
kx- Then the analog of Eq. (6) for the second plate reads 

\lo -vkx\ ^ -uolfexl- (7) 




Figure 1: Friction depends on velocity v through the function g. Below a certain velocity, 



■Wmin = 2t;o = 2c/^/e, the friction force is zero; it starts to rise linearly at Vr, 
maximum and then falls off. 
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The above conditions limit the range of integration to 



kx > 0, 



and 



vokx <uj < {v -vo)kx. 



One then finds the minimum velocity where a frictional force arises as 

c 
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This threshold velocity is reminiscent of the classical Cherenkov effect, although larger by a 
factor of two. However, in the center- of-mass frame where the two plates move at the same 
velocity but in opposite directions, we find the same condition as that of the Cherenkov 
effect: A frictional force arises when, in the center-of-mass frame, the plates' velocity exceeds 
that of light in the medium. As a specific example, we consider a two dimensional space, i.e. 
surfaces represented by straight lines. The dependence of friction force on relative velocity 
is then plotted in Fig. 1. 



3 Formalism and Derivation 

In the previous section we argued for the appearance of friction between moving parallel 
plates which is reminiscent of the Cherenkov effect. Establishing a complete correspondence 
requires a full analysis of the radiation within each object. In this section, we provide several 



arguments to demonstrate why and how a fluctuation-induced friction arises in the context 
of macroscopic objects in relative motion. We start with a heuristic argument, making a 
connection with the Frank- Ginzburg condition. We then present three distinct derivations 
of the friction force using techniques developed in different fields: The first method relies 
on the input-output formalism in a second-quantized picture; the second one appeals to 
quantum field theory in curved space-time. The last, and the longest, derivation is based 
on fiuctuation-dissipation theorem, or the closely related Rytov formalism. The advantage 
of the latter approach is in finding correlation functions inside and outside the two plates 
which can be used to compute the radiation within each plate and in the gap between them. 

3.1 Why is there any friction/radiation? 

To start with, let us consider a space-filling dielectric medium described by a constant real e. 
A wave described by wavevector k satisfies the dispersion relation oj = fo|k|, with vq = cj \fe 
being the speed of light in the medium. This relation describes the spectrum of quantum 
field excitations. If the medium is set in motion, the new spectrum can be deduced simply 
by a Lorentz transformation from the static to the moving frame. Assuming again that the 
speed of light in medium is small (or e is large) , we find 

lli = fo|k| + f A:^; , (10) 

with the medium moving with speed v parallel to the x axis. 

Next consider two (semi-infinite) media one of which moves laterally with velocity f , 
whereas the other is at rest. Although the boundaries modify the dispersion relation, 
we may assume that Eq. (10) approximately describes each medium (with v = for the 
stationary body). This is justified by considering wavepackets away from boundaries. We 
thus have two distinct spectra: The spectrum for the plate at rest is akin to a cone while 
that of the moving plate is tilted towards the positive x axis, as in Fig. 2. For the sake of 
simplicity, we limit ourselves to k = kx^- Let us consider the (spontaneous) production of 
two particles, one in each medium. Since linear momentum is conserved the two particles 
must have opposite momenta (kx and —kx)- This process is energetically favored if the sum 
of the energy of the two particles is negative, that is, spontaneous pair production occurs if 
it lowers the energy of the composite system. This condition is satisfied when 

OJl + UJ2 = 'iVQ\kx\ + vkx < ^ ., (11) 

which is possible only if t; > 2vq . We stress that our argument is not specific to a particular 
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Figure 2: The energy spectra for a medium at rest (solid curve denoted by uji)^ and a 
moving medium (dashed curve denoted by CJ2)) in the {uj^ k^) plane. The spectrum for the 
moving medium is merely tilted. The production of a pair of excitations, indicated by solid 
circles at opposite momenta, is energetically possible for v > 2vo. 

reference frame. If both plates are moving with velocities vi and V2, the velocity v in 
Eq. (11) is replaced by W2 — ^ij thus this equation puts a bound on the relative velocity. 

The above argument is similar to the Landau criterion for obtaining the critical veloc- 
ity of a superfluid flowing past a wall [10]. The instability of the quantum state against 
spontaneous production of elementary excitations (and vortices) breaks the superfluid order 
beyond a certain velocity. Quantum friction provides a close analog to Landau's argument 
in the context of macroscopic bodies. The same line of reasoning is adopted in the work 
of Frank and Ginzburg [4-6] . While this argument correctly predicts the threshold velocity 
for the onset of friction, it does not quantify the magnitude of friction and its dependence 
on system parameters. 



3.2 The input-output formalism 

The input-output formalism deals with the second quantized operators corresponding to in- 
coming and outgoing wave functions, relating them through the classical scattering matrix 
[11-14]. From the (known) distribution of the incoming modes, one can then determine the 
out-flux of the outgoing quanta. The input-output formalism has been used to study the 
dynamical Casimir effect in theoretical [15, 16] as well as experimental [17] contexts, and is 
recently generalized in application to lossy objects by utilizing scattering techniques [18]. 



In the present context of parallel plates, we consider incoming waves as originated well 
within each plate, far away from the gap between them (asymptotic infinity). These waves 
propagate towards the gap, and then scatter (backwards or forwards) to asymptotic infini- 
ties. The incoming and outgoing wave-functions for each plate are normalized such that the 
current density perpendicular to the surface of the plate is unity up to a sign. 

The second-quantized field ^j is now indexed hy i = 1,2 designating the two objects. 
This operator (with implicit index i) is now decomposed into modes 

^ = Vl^ ""''^('^-''ll^-i^ll + '^-k|,C4) + h.c. , (12) 

with two copies of this equation, one for each plate. The wavefunction (pi^]^. , in the object's 
rest frame, is defined as 

1 „ik||-X||itifcx2 



g*K||.X||±*«X^^ W > 0, 



out/in^ r Vfcx _ (13) 

with the prefactor being chosen to ensure the normalization, z measuring the distance from 
the surface, and k± = ^/ew^/c^ — kn. Note that the designation of incoming or outgoing 
for propagating waves depends on the relative signs of oj and k±, as indicated in the above 
equation. For the moving plate, the corresponding wavefunctions are obtained by a Lorentz 
transformation ofoj and kx, while k± (being perpendicular to the velocity) remains invariant. 
With the operators defined above, the input-output relation takes the form 

where S is the 2x2 scattering matrix, and the dependence on uj and kn is implicit. The 
scattering matrix can be straightforwardly computed by matching the wavefunction and its 
first derivatives along the boundaries. Note that a scattering channel relates a wavefunction 
labeled by (c(;,k||) on one plate to (a;',kn), the frequency and wavevector as seen from 
the moving frame, on the second plate. At small velocities, we have a;' ~ w — vkx and 
kn ~ kji. Therefore, a positive-frequency mode on one plate can be coupled to a mode with 
negative frequency on the other plate. However, when the frequency becomes negative, the 
corresponding operator changes character, i.e. an operator a^i^ with negative uj is, in fact, 
a creation operator. We can make this more explicit by writing a^ kn = ol^ _k ^- This 
mixing between positive and negative frequencies is at the heart of the dynamical Casimir 



^The sign of the wavevector is reversed due to complex conjugation. 



effect [18, 19]. In a frequency window where this mixing occurs, i.e. for < a; < vkx, the 
input-output relation is recast as 

aCk|| = -5*11 fliLkii + •S'21 a^^k^uj-k^^ ■ (15) 

One can then compute the expected flux of the outgoing modes, (a°" 'af^^). At zero 
temperature, {a^ af^) = 0, and the only contribution to the outflux is due to the second 
term in the RHS of Eq. (15), resulting in 

(Ck\«rJk„> = 0(^^.-^)i52ip, (16) 

with Q being the Heaviside step function. The radiated energy (per unit time) is then 

V= dwL^-Wdkji hoje{vkx-uj) |52ip. (17) 

Similarly, since friction is the transfer of the lateral momentum, it can be computed by 
replacing Hlo by hkx in the last equation. It is a simple exercise to compute the scattering 
matrix. One can then see that Eq. (17), or, more precisely, its counterpart for friction, is 
indeed in agreement with Eq. (3). The input-output formalism makes the derivation almost 
trivial. 
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Figure 3: a) The operators within each object represent incoming and outgoing modes, 
related in the input-output formalism through the scattering matrix, b) The scattering 
problem is similar to the quantum tunneling over a barrier; friction resulting from transfer 
of momentum by the tunneling quanta. 



Equations (14) and (15) relate annihilation and creation operators which satisfy canon- 
ical commutation relations, 

r^out/in ^out/in ti / \ r /i d\ 

Kc.k|| 'Oi^k,, ']=sgn{u)6ij. (18) 

The function sgn(u;) merely indicates that, for negative frequencies, the creation and anni- 
hilation operators should be identified correctly. The above canonical relations applied to 
Eq. (15) yield 

l-\Su\^ = sgn{oo-vt,)\S2iW (19) 

implying that the scattering amplitude corresponding to the backscattering in the first 
medium is larger than unity for oj < vkx- This is an example of the so-called supers cattering 
due to Zel'dovich [7]: For certain modes, a moving (rotating) object amplifies incoming 
waves indicating that energy is extracted from motion. A closely related phenomenon occurs 
for rotating black holes and is known as the Penrose process [20]. Superradiance signals 
quantum instability of the moving object resulting in energy and momentum radiation, 
and a corresponding exertion of frictional force [7,8,21].^ Equations (17) and (19) can be 
combined to yield 

dwL^-Wdk|, hioe{vkx-oj) (|5iip-i). (20) 

This expression is indeed very similar to quantum radiation from a rotating object (a rotat- 
ing black hole in Ref. [23] or a rotating cylinder in Ref. [21]) with the following substitutions: 
t; — 7- il (linear to angular velocity), kx ^ m (linear to angular momentum), and 5ii to be 
replaced by the scattering matrix of the rotating body. 

It is worth noting that the only contribution to the radiation is from modes with kx > 
oj/v > uj/c, corresponding to evanescent waves in the gap between the two plates. In this 
respect, the radiation is a quantum tunneling process across a barrier (in this case, the gap). 
We can recast the wave equation in a fashion similar to the Schroinger equation as 

(-92 + y)(^ = o, (21) 

with Vi = -(eup' l(? - kn), V2 = -(e('^ - vkxf jf? - k|), and Vgap = \uP' jc? - k||; see 
Fig. 3. The relative motion of the two media results in a steady tunneling of particles of 
opposite momenta from one plate to another, thus leading to the slowdown of the motion. 



An alternative approach to the rotational friction is given in Ref. [22]. 
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3.3 Inner-product method 

To quantize a field theory, a first step is to decompose the quadratic part of the Hamiltonian 
into a collection of (infinite) harmonic oscillators; define the corresponding annihilation and 
creation operators, and impose canonical commutation relations. One can then construct 
the Fock space with the vacuum state of no particles, and excited single- and multi-particle 
states obtained by applying creation operators. In high energy physics, the usual start- 
ing point is empty space, but the above procedure works equally well in the presence of 
background matter, as is the case for the Casimir effect. The reason is that canonical 
quantization only relies upon time translation and time reversal symmetry. The former 
allows construction of eigenmodes of a definite frequency, which is the basis of the notion 
of modes/quanta/particles. Time reversal symmetry, on the other hand, is used to identify 
creation and annihilation operators: The coefficient of a positive (negative) frequency mode 
is understood as an annihilation (creation) operator. To make this correspondence explicit, 
consider a quantum field <]&(t,x), possibly in the presence of a background object which is 
static. One can find a basis of eigenmodes 0aja(x) labeled by frequency co and quantum 
number a to expand the field as 

$(t, x) = yi j; e— Vc.a(x) a^a + e*"V:«(x) oL , (22) 

ci;>0, a 

with (defining 5(x) = 2ttS{x)) 

[au^ia, al^i^l = 5(^1 - L02) 5ai3 ■ (23) 

The latter follows from the canonical commutation relations between the field $(t,x) and 
its conjugate momentum n(t,y), 

[|.(t,x),fl(t,y)]=i;i5(x-y). (24) 

When the object is moving, we lose one or both symmetries in time. The case of 
two parallel plates in lateral motion respects time translation symmetry as the relative 
position does not change. Time reversal symmetry, on the other hand, is broken; in the 
backward direction of time the plate moves in the opposite direction. In the absence of 
time reversal symmetry, the above correspondence between positive (negative) frequency 
and the annihilation (creation) operators breaks down. There is, however, a more general 
way to identify operators as follows. Let us consider two functions (pi and (/)2, which are 
solutions to the classical field equation, and define an inner product as [24] 

(01,02) ='-Jd^i<Pl7T2-7rl<P2), (25) 

11 



where ttj is the corresponding conjugate momentum, and the integral is over the whole space. 
One can easily see that the inner product defined in Eq. (25) is independent of the choice of 
the reference frame or the (space- like) hypersurface as the integration domain. Furthermore, 
we can always find a set of functions (jj^ja which solve the classical field equation and form 
an orthonormal basis, '^ 

{4>uiia, 4>uj2l3) = 5(^1 - ^2)Sal3, (26) 

while their conjugate modes are orthonormal up to a negative sign, 

(Ca,C/3) = -5(^l-^2)<Ja/3. (27) 

These modes form a complete basis, such that the field ^ can be expanded as 

|.(t, x) = yi ^ e-*"Vc.a(x) a^a + e*"V:«(x) aL • (28) 

u)a: positive-norm 

Note that unlike Eq. (22), an annihilation (creation) mode is identified with a positive 
(negative) norm and not frequency; the latter depends on the reference frame while the 
former does not. The above relations become obvious for a field theory in a static back- 
ground where the conjugate momentum is proportional to dt^; Eq. (25) then implies that 
a positive (negative) co corresponds to positive (negative) norm. 

The inner-product method is used in application to quantum field theory in curved 
space-time. For example, it has been employed to show that a rotating black hole is unstable 
due to spontaneous emission [23]. For two parallel plates in motion, we first introduce a 
complete basis. Due to the translation symmetry in time and space (parallel to the plates' 
surface), wavefunctions can be labelled by frequency oj and tangential wavevector kn defined 
in the lab frame where the first plate is at rest and the second one is moving. There are 
two independent solutions defined as 



tJkii 



and 



tjkii 



V^i^k,, + ^11 V^rJkii > plate 1, ^^^^ 

S21^rXn^ plate 2, 



^i2^rJk„, plate 1, ^^^^ 

V'Lk,, + ^22 ^TX„ > plate 2. 



The incoming and outgoing functions are defined in Eq. (13), with the functions in the 
second plate properly Lorentz-transformed; see the explanation below Eq. (13). To find the 



'Note that u> is still a good quantum number because of the translation symmetry in time. 

12 



conjugate momentum, note that Eq. (1) can be schematically derived from a Lagrangian 
C = \\\{dt^f - (V^)^] with VQ = c/v/e, hence H = dL/d^ = \dt^. Similarly for a 
moving object C ss ^[^(dt^ + ydx^)"^ - (V<1>)^] and 11 = -^{dt^ + vdx^)- In terms of 
partial waves, -/r^jk = —i ^~'"i "" 0^k within the moving plate and similarly for the static plate 



with V = 0.^ One can then see that 

{^l,^ , 'AL.i,, ) = 3(^1 - ^2) 5(k|| - 1|| ), ^1 > 0, (31) 

and 

('/'"k||,'/'"i||) = sgn{uJi-vk^)6{coi-co2)3m - 1||), u;i > 0. (32) 

(We have exploited the fact that the norms are diagonal in frequency to obtain the delta 
functions in k± which are then converted to those of uj.^) Note that the (super) unitary 
relation in Eq. (19) is essential in deriving the norms. Functions of type I have positive 
norm so they serve as the coefficients of annihilation operators. However, type II functions 
include negative- norm modes for < a; < vkx- Therefore, despite the positive sign of 
frequency, the latter should be identified as creation operators. We thus expand the field as 
l>(t,x) ^ 

i^>0,kii 



^^/2 E e"^"VLk|,(x)aLk„ + 



J2 e— Vl,|(x)al,|+ E e— Vi}k|,(x)aL^kt +h.c.,(33) 

0<a;,tifca;<w ,k|| 0<w<t!fca;,k|| 

where the summation is a shorthand for multiple integrals, and a and a' satisfy the usual 
commutation relations. The radiation within the first plate is given by 

= |^^<lu|d^^{-{l-\Sn\') + \S^,\'} 

f'OO f' 

= acj dk\ihuje{vkx-u;)\Si2\'^, (34) 

where we have again exploited Eq. (19), arriving at the same results as in the previous 
sections. Notice that only superradiating modes contribute to the radiation while other 
modes cancel out in the second line of the last equation. 

3.4 Radiated energy: The Rytov formalism 

In this section, we employ the Rytov formalism [9] to study the correspondence between 
friction and radiation in some detail. This formalism is based on the fluctuation-dissipation 



^An alternative is to define the integration hypersurface to consist of cacli object's rest frame. 

®The integral over the gap does not contribute to the frequency delta functions and thus is neglected. 
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theorem, and has been extensively used in the context of non-contact friction [2]; see also 
Ref. [3] and references therein. This section goes beyond the existing literature by computing 
various correlation functions and the energy radiation within the plates (as well as in the 
gap between them) and discuss their dependence on the reference frame. 

We start by relating fluctuations of the field to those of "sources" within each medium 
by 

A + ^e(w,x)j$(w,x) = -^/9^(x). (35) 

The "charge" p fluctuates around zero mean with correlations (co- variance) 

{Pu:{^)plj{y)) = a(w)Ime(tj,x)(5(x-y), (36) 

where 



a{uj) = 2h 



"(w,T) + - 



.coth(J^), (37) 



The field is related to the sources via the Green's function, G, defined by 

A + !^e(w,x)J G(w,x,z) = (5(x-z). (38) 

In equilibrium (uniform temperature and static), this results in the field correlations 
(a>(a;,x)$*(u;,y)) = ^ / dz G(o;,x,z)G*(o;,y, z) (/,^(z)p:(z)) 

'-- ^AU snace 



a;2 



2 a(ti;) / dzG(a;,x, z) Ime(aj, z) G*(ti;,y, z) 

J All space 



= a(u;) ImG(a;,x, y), (39) 

in agreement with the fluctuation-dissipation condition which relates correlation functions 
to dissipation through the imaginary part of the response function. Note that the second 
line in Eq. (39) follows from ^Ime = — ImG~^ according to Eq. (38). 

For the case of two plates, we first compute the correlation function for two points in the 
gap. The source fluctuations in each plate will be treated separately; starting with those in 
the static plate (indicated by subindex 1 on the integral) 

(<5(w,x)<5*(w,y))i = -^ai{uj) / dzG{io,^,z)lme{uj,z)G*{uj,y,z) 
c Ji 

= ^^ai(^) / dz [e(a;,z)G(tj,x,z)] G*(w,y,z) - G(u;,x,z) [e{uj,z)G{uj,y,z)]* 

= -ai(a;) / dz [AzG(a;,x,z)] G*(a;,y,z) - G(a;,x,z) AzG*(a;,y,z) 

= ^ai(w) /dS-[(VzG(^,x,z))G*(^,y,z)-G(L^,x,z)V,G*(^,y,z)] . (40) 
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Note that we used Eq. (38) in going from the second to the third hne above, and then inte- 
grated by parts to obtain an integral over the surface adjacent to the gap. The contribution 
due to the other surface at infinity vanishes since e is assumed to have a small but finite 
imaginary part. 

To compute the surface integral in Eq. (40), one needs the (out-out) Green's function 
with both points in the gap. The latter is given by Eq. (A.l), and leads to 

„ (. ,\ \ ipad\2 

($(L^,x)cI>*(^,y))l = -y^iM '"^ ' , \U{Ra)\''x 



In this equation, R^ (Ra) is the reflection coefficient from the first (second) object with a 
being a shorthand for both the frequency and wavevector. Also x/y (x/y) is the distance 
from a reference point on the surface of the first (second) plate — the reference points on two 
surfaces have identical parallel components, xn = xii , but differ in their z-component as z^ = 
d — Zy^. The regular and outgoing functions are defined with respect to the corresponding 
plate; see Appendix A for more details. Furthermore, the overbar notation implies complex 
conjugation, and \U{Ra)\^ is defined as 

\JJ(J^ )|2 = JdS-[$^V$.-(V<^*)c|>.] 



- Q 



with $Q, = $a^(a;, z) + i?Q,$°"*(a;, z), such that 

\U{Ra)\ = 1 — \Ra\ , d G propagating waves, 

= 21m.Ra, a G evanescent waves. (43) 

One can similarly find the correlation function due to source fiuctuations in the second plate 

(<^( uj, x)<^* ( oj, y))2 = -y -^ — : — — \U(Ra)\x 



(C^(a;,x) + i?a$r(^,x))(C'(a;,y)+i?^,<I>r(^,y)). (44) 

The total correlation function is the sum of the contributions due to each plate, 

($((^,x)$*(w,y)) = (cD(^,x)cI>*(L^,y))i + ($(l^,x)$*(c^, y))2. 
The frictional force is then computed as the average of the appropriate component of the 
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stress tensor as 



/OO f' 

OO J 

/* OO /* \ 7 fc \ (i \ Zi 

= / duL""-' d^hk J" I |[/(i^^i,„)|^|^(i;^i,„)|=^(niM-n2(a;-vk)), 

JO J |i e ^ -K(^k|| -K(^k|| I 

(45) 

where we have restored kn in place of a. Further manipulations lead to Eq. (3). Similarly, 
the Rytov formalism allows us to calculate the energy flux from one object to the other as 

/OO f 

doj / dS(5i$(w,x)a,$*(^,x)) 
OO J 

/* OO /* I 7 A* I /Y 1 2 

= / dioL^-^d^ fu^ J^\ J |^(i?^i,„)p|[/(fi^i,„)p(niH-n2(a.-vk)), 

(46) 

i.e. by merely replacing hkx with fwj in Eq. (45). 

Next we compute the energy flux through each plate. Since the dielectric function is 
assumed to be a real constant (albeit with a vanishingly small imaginary part), we can 
circumvent ambiguities in defining the Maxwell stress tensor in a lossy medium [25]. In 
the following, we find the field correlation function in the first (static) plate due to source 
fluctuations in the moving plate, by using an analog of Eq. (40) but evaluating a surface 
integral on the second plate. However, in this case, the appropriate Green's function is the 
(out-in) type given in Eq. (A. 4). We then find the latter correlation function as 

($(c.,x)$*(.,y)), = y -^i----^'^) k^^T^ i^^^^p^ 



(y„Va(w,x) +H^„Va(c^,x)) (14^^(6^, y) +VFaVa(c^,y)) , (47) 

where x and y are both inside the first plate, V and W are coefficients depending on a and 
system parameters, and the functions il) are defined inside the medium; see Appendix A for 
more details. Henceforth, we assume that the objects are at zero temperature. Anticipating 
that only evanescent waves contribute, we obtain the energy flux in the first plate due to 
the fluctuations in the second plate as ^ 

Pi^) = /d.;L^-i/"dk,| hu: — — ;g|^^^^_- 2Imi?.k,, 2Imi?.k„ sgn(a; - v k), 

(48) 
^Inside the object, the "Poynting vector" is defined by the same tensor dt'^dz'^. 
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where we have used the fact that, for evanescent waves, Pa = k± = ./co'^/c'^ — kn is purely 
imaginary while pa = k± = ^/ew^/c^ — kn is real, leading to 



\VJ'^-\WJ^ = —21mRa. (49) 

Pa 

In order to take into account the source fluctuations in the first plate (where we compute the 
field correlation function), we need the (in-in) Green's function in Eq. (A. 7). The energy 
flux due to the latter fluctuations, P^ , is computed similarly but there is one subtlety. 
Unlike the previous cases, the correlation function is evaluated at points where there are 
also fluctuating sources. However, Eq. (40) contains, beyond the surface integral, a term 
proportional to ImG((x',x, y) which does not contribute to the radiation. The remaining 
computation is similar to the previous case, and the overall energy flux is obtained as 

p(i)=p«+p« 

doj L^-7 dk|, huj ^^^^ ^ e(v k - u). (50) 

This is again in harmony with the results in the previous sections. 

Comparing Eqs. (50) and (46), we observe that in the reference frame in which the first 
plate is at rest, 

P^'^ = ^gap- (51) 

However, Pgap must vanish in the center of mass (cm.) frame from symmetry considera- 
tions. It can be obtained explicitly by a Lorentz transformation from the lab frame, which, 
to the lowest order in velocity, takes the form 

= P|-- = Pgap-^//2, (52) 

indicating P^^' = vf /2. This conclusion can be verified directly as follows: First note that 
Eqs. (45) and (50) yield 

pa._^./ ^,„r%.4._^v-*'^^'-'^"-""'^\ (53) 

2 A,>0 'Jo \ 2 J |l-e2*-f'i;„k,«,t,|2 ' 

where the x axis is chosen parallel to the velocity v. Let us make the following change of 
variables 

Lo' = u — vkx/2, 

k'x = kx — vuj/2c ss kx, 

k'i = ki, i y^ X. (54) 
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It then follows that 

p(i) --^= dk / dL^'W — ^, ' '' — ^ , (55) 

^ ifei>0 J-vk'^/2 |i e i ^a;',k]|-"a;',k||l 

where i?"'' and i?~ are the reflection matrices from plates moving at velocities v/2 and —v/2, 
respectively, along the x axis. Since e is real^, we have 

R+ ,,, = R-, . , . (56) 

This implies that the integrand in Eq. (55) is antisymmetric with respect to io' so that the 
integral vanishes. 

When there is friction, work must be done to keep the moving plate in steady motion. 
This work should be equal to the total energy dissipated in the plates, 

vf = Ptot, (57) 

where Ptot is the sum of energy flux through each plate. For Eq. (57) to hold, the energy 
flux through the second (moving) plate should also be equal to P^^' = vf /2. In the center- 
of-mass frame too, we should have the same condition because of the energy conservation 
vf = Pc.m. +Pc.m., and the symmetry PCm. = -fc.m.- (The force in the center- of -mass frame 
is almost identical to the lab frame since the velocity is small compared to the speed of 
light.) Therefore, we conjecture that P^ = Pg = vf/2 irrespective of the reference frame 
S, while PgUp sensitively depends on the reference frame S; it is vf/2 when the first plate 
is at rest, —vf/2 when the second plate is at rest, and zero in the center of mass frame. 



4 Discussion and Summary 

Throughout this paper, we explicitly considered semi-infinite plates described by a constant 
and real dielectric function. However, the underlying physics is rather general, and does not 
depend on the idealizations that made for the sake of convenience. For example, rather than 
a semi-infinite plate, we can consider a thick slab of a material with a complex dielectric 
function e. The slab will act like an infinite medium provided that the imaginary part of e 
while small, is sufficiently large to absorb the emitted energy within the slab, with almost 
no radiation escaping the far end. Such conditions can be met for a broad range of the 
thickness and lossy- ness. 



In general, the real part of the response function is even in frequency, i.e. Ree(cj) — Ree(— oj) 
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Classically, Cherenkov radiation is emitted when a charged particle passes through a 
medium. However, even a source without a net charge, or even a multipole moment, may 
result in Cherenkov radiation due to quantum fluctuations [4]. In the present paper, this is 
demonstrated for two neutral parallel plates in relative motion. By employing an amalgam 
of techniques, usually applied in different contexts, we are able to make novel conceptual 
and technical observations. These techniques are applicable to a variety of other setups: 
An interesting situation, closer in spirit to Cherenkov radiation, is when a particle passes 
through a small channel drilled into a dielectric. Another closely related problem is a 
particle moving parallel to a surface [18, 26, 27]. Our approach of utilizing scattering theory 
in conjunction with a host of other methods, including input-output and Rytov formalism, 
should be useful in analyzing such situations. 
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A Green's functions 

In this appendix we compute a number of Green's functions where the two spatial arguments 
lie within or outside each plate. We take the flrst plate to be at rest while the second one is 
moving at a velocity v parallel to its surface. We further assume that |v| ^ c for simplicity. 
• Green's function with both points lying within the gap (outside both objects): In this case, 
the Green's function is given by (with z^ > z^,) 

1 p^P°''^ 

Gout-out(a;,x,z) = j;--:^ (.I>gg(.;,x)+fi^.^r(^,x))x 

($-s(a;,z) + i?„$r(^,z)), (A.l) 
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where we used a compact notation defined as 



a = k|| , a = — k|| , Pa = k± = Joj'^ — kn , 

^reg _ ik||-X||+ifcx2 ^out _ Jk||-X||-iA:x2 

^a ^ 1 ^a ^ 1 

Ra = -Ri(a;,k||) = i?^k|| , Ra = i?2(^,k||) = i^^/kj^ , 

zx = rf-^x, x||=x||, ^ = L^-iy"dk||^L^-i /^-^, (A.2) 



iD-U 

(2T) 



where w' and kj, are the Lorentz transformation of uj and kn, respectively. Also D is the 
number of (spatial) dimensions. According to these definitions, <&a(z) is defined with respect 
to an origin on the surface of the first plate, while <^a(x) is the wavefunction defined with 
the origin on the surface of the second plate and the direction of the z-axis reversed. It is 
straightforward to check that the expression in Eq. (A.l) is indeed the Green's function. 
First note that, for x 7^ z, it solves the homogenous version of Eq. (38). Furthermore, the 
coefficients are chosen to produce a delta function when z — )■ x upon applying the Helmholtz 
operator. 

• Green's function with one point in the gap and the other inside the first plate: This 
Green's function can be obtained from continuity conditions, i.e. by matching the Green's 
functions approaching a point on the boundary from inside and outside the object 

G'out-out(w,x,y)|y^5.+ = Gout-in(^,x,y)|y^2-- (^-3) 

This leads to 

piPad 

Gout-in(u;,x,z) = ^l/(2ip^) (cDg^g(a;,x) + i?^c&r(a;,x))x 

Q, i e ™ itctiXa 

(y«V+(^,z) + T^aC(^,z)), (A.4) 

with 

^±^gik||.x|i±ifc^2^ (A.5) 



where k± = pa = A/ew^/c^ — kn. The (diagonal) matrices V and W are determined by 
imposing continuity equations, as 

Va + Wa = l + Ra, 

PaiVa - Wa) = Pail " Ra) ■ (A.6) 
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• Green's function with both points inside the first plate: This Green's function is given by 

Gin-in (W, X, Z) = 2^ l/(2ipa) ^ (KvV'a(^. X) + TV^^" (oj, x)) X 

(T4^t(a;,z) + I^aV^(w,z)), (A.7) 

where, via continuity relations, we have 

PaiVa - Wa) = Pc.e'^^-\l " Ra)- (A.8) 
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